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Circle  Detection  for  Extracting  Eddy  Size  and 
Position  from  Satellite  Imagery  of  the  Ocean 

Sarah  H.  Peckinpaugh  and  Ronald  J.  Holyer 


Abstract— This  paper  evaluates  the  ability  of  several  circle  de¬ 
tectors  to  define  the  size  and  position  of  warm  and  cold  eddies 
in  oceanographic  satellite  imagery.  The  Advanced  Very  High 
Resolution  Radiometer  channel  4  (10.3-11 .3  pm)  or  sea  surface 
temperature  test  images  are  reduced  to  binary  edge  images.  Six 
different  circle  detectors  are  then  applied  to  the  edge  images. 
The  automated  results  are  compared  to  eddies  defined  by  a 
trained  analyst. 

Introduction 

NFRARED  (IR)  images  of  the  ocean  obtained  from  sat¬ 
ellite  sensors  are  widely  used  for  the  study  of  ocean 
dynamics  [  1  ]— [6] .  Many  studies  involve  analyzing  large 
data  sets  of  IR  imagery.  It  is  desirable  to  replace  the  la¬ 
bor-intensive  time-consuming  manual  interpretation  of  IR 
imagery  with  automated  analysis.  And,  because  of  the 
current  worldwide  awareness  of  the  effects  of  climatology 
and  global  change,  the  impetus  is  now  on  both  speed  and 
accuracy. 

Therefore,  the  Naval  Research  Laboratory  began  de¬ 
velopment  of  the  Semi-Automated  Mesoscale  Analysis 
System  (SAMAS)— a  comprehensive  set  of  algorithms  that 
handles  the  entire  automated  analysis  problem.  The  SA¬ 
MAS  performs  the  low-level  edge  detection  through  fea¬ 
ture  formation  and  higher  level  artificial  intelligence  mod¬ 
ules  that  estimate  positions  of  previously  detected  features 
when  cloud  cover  obscures  direct  observation  in  the  cur¬ 
rent  image  set  [7],  Feature  locations  produced  automati¬ 
cally  by  SAMAS  compare  favorably  to  human  interpre¬ 
tation  of  the  same  images  [8],  [9].  An  important  part  of 
the  SAMAS  is  the  eddy  detection  based  on  an  ensemble 
of  edge  fragments  that  have  been  labeled  by  a  feature- 
labeling  module  as  being  part  of  an  eddy  [10].  The  SA¬ 
MAS  uses  the  Hough  Transform  to  define  eddy  position 
and  size.  This  paper  evaluates  the  effectiveness  and  com¬ 
pares  the  results  for  six  methods  of  circle  detection,  in¬ 
cluding  the  Hough  Transform  method. 

Preprocessing 

The  test  data  consist  of  1 8  IR  images  of  the  Gulf  Stream 
area  of  the  North  Atlantic  Ocean.  These  images  were  ob- 
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tained  from  the  AVHRR  sensor  aboard  a  National  Oceanic 
and  Atmospheric  Administration  satellite.  Images  are  cal¬ 
ibrated;  the  channel  4  images  are  scaled  to  8  bits,  such 
that  0  corresponds  to0°C  and  255  corresponds  to  25.5°C, 
and  are  then  remapped  to  a  Mercator  projection.  Each  im¬ 
age  is  256  x  256  pixels  with  a  warm  or  cold  eddy  in  the 
center  of  the  image.  Resolution  of  the  test  images  is  about 
2.5  km/pixel.  From  each  image,  a  binary  edge  image  was 
created. 

The  method  of  edge  detection  used  to  create  the  edges 
is  described  by  Holyer  and  Peckinpaugh  [II],  but  we 
modified  the  zero  crossing  test  slightly.  The  modification 
uses  a  second  threshold  to  extend  edges  found  using  the 
original  method.  The  results  of  this  type  of  edge  detection 
are  2-pixel-wide  lines  at  the  edges.  The  lines  are  dilated 
from  2-pixel-wide  edges  to  4-pixel-wide  edges.  An  ana¬ 
lyst  subjectively  edited  the  edges  and  kept  only  those  con¬ 
sidered  to  be  associated  with  the  eddy.  A  second  analyst 
looked  at  only  the  IR  images  using  an  interactive  editor, 
and  subjectively  defined  the  center  coordinates  and  the 
radius  of  the  eddies  for  each  image.  Fig.  1  shows  the  test 
images,  the  detected  edges,  and  the  circle  fit  subjectively 
defined  by  the  human  expert. 

These  visually  defined  center  coordinates  and  radii  will 
be  used  to  evaluate  the  eddy  detection  methods.  In  Fig. 
1,  darker  shades  indicate  warmer  temperatures,  red  indi¬ 
cates  edges,  and  yellow  indicates  the  analyst-defined  ed¬ 
dies.  The  eddy  images  selected  as  test  cases  are  typical  of 
those  encountered  in  a  day-to-day  operational  environ¬ 
ment.  That  is,  some  have  very  weak  surface  expression 
and  others  are  partly  cloud  covered.  The  test  has  not  been 
biased  by  inclusion  of  only  ideal  test  cases.  Cold  eddies 
are  generally  more  difficult  to  observe  than  warm  eddies 
because  the  dense,  cold  water  tends  to  sink  quickly, 
thereby  greatly  reducing  the  surface  IR  signature.  The  test 
cases  include  10  cold  eddies  [Fig.  1  (a)— ( j)]  and  8  warm 
eddies  [Fig.  l(k)-(r)]. 


Circle  Detection 

The  six  algorithms  selected  for  testing  are  referred  to 
by  the  following  names;  Hough,  Thomas_Chan,  Landau, 
Albano,  Janowitz,  and  Circle  Fit.  Some  of  these  algo¬ 
rithms  were  modified  slightly  from  the  original  descrip¬ 
tions  to  tailor  the  algorithm  better  to  the  input  data  or  to 
the  desired  output.  The  methods  as  they  were  used  in  this 
study  are  defined  briefly  below. 
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Fig.  I.  IR  images  of  eddies.  Actual  edges  are  red.  and  analysis-defined 
eddies  are  yellow,  (a)  Cold  eddy,  SST  image  centered  at  36.09  latitude 
-71.18  longitude;  Jan.  14,  1987.  (b)  Cold  eddy,  SST  image  centered  at 
35.89  latitude  -70.53  longitude;  Mar.  7.  1987.  (c)  Cold  eddy,  SST  image 
centered  at  36.46  latitude  -67.24  longitude;  May  11,  1987.  (d)  Cold  eddy. 
SST  image  centered  at  34.95  latitude  -66.29  longitude;  June  1.  1987.  (e) 
Cold  eddy.  SST  image  centered  at  35.69  longitude  -68.34  longitude;  June 
26,  1986.  (f)  Cold  eddy,  SST  image  centered  at  35.89  latitude  -67.94 
longitude;  July  8,  1986.  (g)  Cold  eddy.  SST  image  centered  at  37.82  lati¬ 
tude  -66.69  longitude:  July  15,  1986.  (h>  Cold  eddy,  SST  image  centered 
at  33.75  latitude  -71.68  longitude:  May  11,  1987.  (i)  Cold  eddy,  SST 
image  centered  at  36.18  latitude  -63.40  longitude;  June  11.  1987.  ( j) 
Cold  eddy.  SST  image  centered  at  37.62  latitude  -57.66  longitude:  May 
11.  1987.  (k)  Warm  eddy,  Channel  4  image  centered  at  40.23  latitude 
-63.55  longitude:  Dec.  4,  1986.  (1)  Warm  eddy,  SST  image  centered  at 
38.88  latitude  -70.13  longitude;  Mar.  7.  1987.  (m)  Warm  eddy.  SST  im¬ 
age  centered  at  40.80  latitude  -57.31  longitude;  May  1 1,  1987.  (n)  Warm 
eddy,  SST  image  centered  at  41.74  latitude  -58.31  longitude;  June  II. 
1987.  (o)  Warm  eddy,  SST  image  centered  at  39.46  latitude  -68.04  lon¬ 
gitude;  Sept.  5,  1986.  (p)  Warm  eddy,  SST  image  centered  at  39.34  lati¬ 
tude  -70.03  longitude;  Nov.  17,  1986.  (q)  Warm  eddy,  Channel  4  image 
centered  at  39.93  latitude  -65.54  longitude;  Mar.  25.  1987.  (r)  Warm 
eddy,  SST  image  centered  at  39.07  latitude  -68.04  longitude;  Apr.  11, 
1987. 


Hough 

Duda  and  Hart  [12]  provide  a  general  description  of  the 
circular  Hough  Transform.  The  Hough  method  can  be 
conceptualized  as  a  convolution  of  various  binary  kernels, 
each  containing  a  circle  of  “  Is”  of  some  selected  radius, 
with  the  binary  edge  image.  If  all  possible  radii  are  rep¬ 
resented  in  the  library  of  convolution  kernels,  the  kernel 
resulting  in  the  highest  convolution  result  identifies  the 
best  estimate  of  radius  for  any  circular  structure  that  might 
be  in  the  edge  image.  The  x  and  y  convolution  lag  values, 
where  the  maximum  output  occurs,  identify  the  best  es¬ 
timate  of  center  coordinates  for  the  circular  structure.  This 
convolution  procedure  can  be  shown  to  be  mathematically 


optimal  in  the  least-mean-squares  sense  for  circles  con¬ 
tained  as  additive  signal  in  a  random,  uncorrelated  noise 
background.  The  mathematical  optimality  may  not  strictly 
apply  in  the  present  case  where  circles  are  represented  by 
scattered  edge  fragments.  However,  the  detection  per¬ 
formance  known  to  result  from  convolution  of  an  image 
with  a  kernel  containing  the  desired  target  suggests  that 
the  Hough  Transform  method  should  be  investigated. 

In  practice,  the  Hough  method  does  not  perform  the 
multiple  convolutions  described.  Rather,  it  performs  an 
equivalent  but  much  more  efficient  operation,  which  in¬ 
volves  incrementing  the  elements  of  an  accumulator  array 
for  each  edge  pixel  in  the  edge  image.  The  final  accu¬ 
mulator  array  is  then  searched  for  the  maximum  value. 
The  address  of  the  maximum  value  in  the  accumulator 

•ay  defines  center  coordinates  and  the  radius  of  a  circle, 

...resenting  a  best  fit  to  the  image  edge  fragments. 

,7Tr  the  sake  of  better  computational  efficiency,  our  im- 
t  .ementation  of  the  Hough  method  uses  a  two-dimen¬ 
sional  accumulator  rather  than  a  three-dimensional  array. 
The  third  dimension  is  avoided  by  not  simultaneously 
considering  all  possible  radius  values,  but  by  considering 
one  radius  value  at  a  time.  The  starting  point  for  radius  is 
the  largest  geophysically  reasonable  radius  (133  km  for 
this  study).  Radius  is  then  decremented  until  a  case  is 
reached  where  a  circle  is  found  widch  has  40%  of  its  cir¬ 
cumference  occupied  by  an  edge  pixel. 

Given  a  set  of  points  defined  as  (xh  y, ),  /  =  1 , 2,  •  •  •  , 
N,  the  Hough  Transform  is  applied  as  follows. 

a)  Convert  the  upper  limit  of  1 33  km  to  a  discrete  pixel 
measure,  R, . 

b)  Zero  all  values  of  the  two-dimensional  accumulator 
array,  A. 

c)  Using  Bresenham's  [13]  algorithm,  compute  the  list 
of  discrete  coordinates  ( xh ,  yh)  for  a  circle  of  radius 
Rj  centered  at  (0,  0).  Also  computed  is  the  number 
of  points  forming  the  circle,  Nh. 

d)  Increment  values  of  A:  Ax,  +  ,„.v, ,.+  =  Ax,  +  n 

+  1 . 

e)  Select  the  maximum  value  from  A ,  Axm  ym.  If  ( Nh  * 
0.4)  >  Axm  ,.m,  then  RJ  + ,  =  /?,  -  1,  and  go  to  step 
b. 

f)  Xc  =  xm,  Yc  -  ym,  and  R  =  Rj. 

Thomas  jChan 

Given  a  set  of  points  defined  as  (jc;,  y,  ),  i  =  1,  2,  •  •  •  , 
N,  Thomas  and  Chan  [14]  propose  that  the  center  position 
and  radius 

x  =  C|h2  -  c2b]  =  Q| e2  -  a2f| 

a\b2  —  a2b\  c  a{b2  -  a2h, 

can  be  computed  as  follows: 

-  2  I \y,r,  +  Nr-  'j 
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where 


at  =  2\  Yxj  —  NYx2 


b{  -  2  Xi  Y  y,  -  N  Y  Xi  y^j 
a2  =  2  (Yi  Xj  Y  yt  -  N  Y  xt  y, \  =  b 


c,  =  ^S  x2  Y  x, ■  -  N  Y  xf  +  Y  Xi  Y  yj 
-  N  Y  Xjyf'j 

c2  =  ( Y  x]  Y  y,  -  N  Y  yf  +  Y  y,  Y  y2 


-  N  YxJy^J 


These  equations  result  from  minimization  of  a  theoretical 
error  function  in  Xc,  Yc,  and  R.  The  reader  is  referred  to 
the  original  paper  for  the  definition  of  the  area-based  error 
function  and  the  mathematics  of  minimization  leading  to 
the  equations  listed  here. 


Landau 

Given  a  set  of  points  defined  as  (*,,  y,),  /  =  1,  2,  •  •  •  , 
N,  Landau  [15]  suggests  an  iterative  method  for  comput¬ 
ing  the  center  position  and  radius.  The  algorithm  is  ap¬ 
plied  to  the  data  as  follows. 

a)  For  a  =  (Xc,  Yc),  the  center  coordinates  of  the  cir¬ 
cle,  initialize  aQ  =  {Xa,  Ya),  where 

Xa  '  N  ^  Ya  ~  N  ^  y‘- 

b)  Use  5*  to  compute  the  associated  radius  Rk 

Rk  =  ^Y\?i  -  dk\ 
where  r,  =  (x„  y(). 

c)  Use  a*  and  Rk  to  compute  ak  + 1 

ak+ 1  =  Y  (?,•  —  Rkin~at). 

d)  If  |  a  —  a*  + 1 1  >  0.000001,  then  go  to  step  b. 

This  method  also  involves  minimization  of  an  error 
function,  which  is  the  LMS  distance  from  the  set  of  data 
points  to  the  circle  defined  by  Xc,  Yc,  and  R.  Minimization 
of  the  LMS  error  function  does  not  lead  to  an  exact  so¬ 
lution,  thus  the  need  for  the  iterative  approach. 


Albano 

Albano  [16]  suggests  solving  the  following  general 
conic  equation: 

Ax2  +  Bxy  +  Cy 2  +  Dx  +  Ey  +  1.0  =  0 

to  fit  conic  sections  to  edge  points. 

This  method  is  more  general  than  the  others  in  that  it 
can  be  used  to  fit  any  conic  arc,  not  just  circles,  to  the 
data.  However,  we  have  constrained  the  solution  here  to 
provide  only  circular  fits.  The  method  employs  least- 
squares  minimization  of  the  conic  equation  evaluated  at 
ail  data  points.  Coordinate  standardization  [16]  is  also 
performed  in  conjunction  with  this  method. 

The  steps  to  compute  the  center  coordinates  ( Xc ,  Yc) 
and  radius  R  are  as  follows. 

a)  Coordinate  standardization: 

Xb  =  N  Y  Xi  Yh  =  N  2  y‘ 

x,  =  (x,  -  Xb)  y,  =  ( >-,  -  Yh) 

S  =  Vs  Xj  +  Y  y, 


b)  Solve  the  following  simultaneous  linear  equations 
for  A,  B,  C,  D,  and  E. 

A  S-4  +  B Yxl ysi  +  CY x 2  y2  +  DYxl 

+  E  Yx2,  ysi  +  Yxl  =  0 

A  Yx3si ysi  +  BYxl yl  +  C Yxsl  yl  +  D  Yx2,  ysi 
+  EYxsiyl  +  Yxsiysi  =  0 
A  Yxl  yl  +  B  Yxsi  yl  +  CYyl 

+  D  Y xsi  y2si  +  E  Y yl  +  S yl  =  0 
A  Yxl  +  B  Yxl  y*i  +  C  Yxsl  yl 
+  D  Yxl  +  E  Yxsi  ysi  +  Yysi  =  0 
A  Y  x2ysi  +  B  Yxsiy2si  +  C  S  yl 
+  D  Yxsiysi  +  E  Yy2  +  Yys,  =  0. 


c)  Compute  the  center  coordinates: 


2  CD  -  BE 
B 2  -  4 AC 


S  +  Xh 


(2AE  -  BD\ 

Yc  \B2-4AC)S  Yb' 

d)  Compute  the  radius  for  a  circle: 

_  £  J(Xr  -x,)2  +  (Yc~  y,)2 
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Janowitz 

The  Janowitz  method  has  been  modified  from  an  ellipse 
detection  scheme  [4]  to  a  circle  detection  scheme.  Jan¬ 
owitz  also  did  a  proximity  test  of  edge  pixels  to  a  candi¬ 
date  center  pixel  position.  If  an  edge  pixel  was  deter¬ 
mined  to  be  close  to  the  candidate  center  pixel  position, 
then  that  position  was  no  longer  considered  a  candidate. 
This  test  did  not  work  well  with  these  edges  and  was  omit¬ 
ted. 

This  algorithm  shoots  rays  in  52  directions  from  each 
candidate  center  point  and  records  the  distance  along  the 
ray  to  the  first  encounter  of  an  edge  pixel.  The  sum  of  the 
distances  between  the  candidate  circle  and  these  ray  dis¬ 
tances  is  calculated,  and  the  minimum  value  found  indi¬ 
cates  the  best  fit.  Since  the  method  considers  only  the  first 
edge  point  encountered  along  each  ray,  one  would  expect 
the  Janowitz  method  to  be  biased  toward  smaller  radius 
values. 

Given  a  binary  edge  image  defined  as  f(x,  y),  the  center 
coordinates  (X0  yc)  and  radius  R  are  computed  as  follows. 

a)  Apply  edge  thinning  algorithm  to  the  image,  / ( x ,  y) 
[17]. 

b)  Initialize  variables 

*0  =  Xc  =  *cm 

Yo  =  Yc  =  Ycm 

MinErr  =  9999.0 

where  xcm,  ycm  are  the  first  possible  coordinates  of 
the  image  /  (jc,  y )  that  can  be  considered  as  a  pos¬ 
sible  center  coordinate. 

c)  Given  the  set  of  52  evenly  incremented  rays  ema¬ 
nating  from  the  center  image  coordinates  X0,  Y0  to 
the  outer  edge  of  the  window  centered  about  that 
point,  find  the  distances  dt  (j  =  1,  2,  •  •  •  ,  52) 
from  the  window  center  to  the  first  edge  pixel,  if  it 
exists,  located  along  the  ray. 

d)  MinRad  =  MIN(33,  dt,  d2,  •  •  •  ,  d52 )  and  MaxRad 
=  MAX(33,  du  d2,  •  •  •  ,  d52).  Compute  the  error 
statistic  for  each  radius  between  the  limits 

S  \(dj  -  m) | 


to  =  MinRad,  •  •  •  ,  MaxRad, 

where  only  those  dj  values  are  used  for  which  an 
edge  pixel  was  encountered  along  the  ray  and  within 
the  window,  and  c  is  the  count  of  those  dj  values 
actually  used  to  compute  the  statistic.  If  (c  <  (0.4 
*  52)),  then  Em  =  9999.0. 

e)  If  (MIN(£MinRad,  •  •  •  ,  £M«xRad)  <  MinErr),  then 
MinErr  =  MIN(£MinRad,  •  •  •  ,  £Ma)(Rad) 

Xc  =  Xo 
Ye  =  Y0 

R  =  radius  associated  with  the  new  MinErr. 


f)  If  Xq,  Y0  are  not  the  last  possible  center  coordinates 
within  the  image  / (x,  y),  then  get  next  X0,  T0,  and 
go  to  step  c. 

CircleFit 

CircleFit  used  the  mean  (centroid)  of  the  edge  points 
as  the  circle  center,  and  the  mean  distance  of  the  points 
from  the  center  as  the  radius.  Given  a  set  of  points  defined 
(xh  y,),  i  =  1,  2,  •  •  •  ,  N,  the  center  coordinates  (X,.,  Yt) 
and  the  radius  R  are  computed  as  follows: 

„  2  V(X(.  -  JC,)2  +  (Yc  -  y,)2 


Experimental  Results 

Figs.  2  and  3  are  the  results  of  the  six  methods  applied 
to  the  18  test  images.  Some  methods  produced  no  valid 
definition  for  some  of  the  eddies.  These  no-definition 
cases  are  not  used  for  the  analysis  of  that  method.  The 
results  of  the  automated  methods  are  compared  to  the  an¬ 
alyst’s  definition  of  the  eddy. 

For  each  eddy,  errors  in  center  position  and  radius,  rel¬ 
ative  to  human  judgment,  were  calculated  for  each  of  the 
circle-detection  methods.  The  circle-detection  methods 
were  then  compared,  pairwise,  based  on  the  mean  posi¬ 
tion  and  radius  errors.  The  statistical  validity  of  the  com¬ 
parison  was  established  using  a  T-test  [18].  The  T-test 
used  18  samples  (or  less,  in  the  cases  having  no-definition 
for  an  eddy)  and  two  treatments.  The  results  of  the  com¬ 
parisons  are  shown  in  Tables  I  and  II,  where  the  best 
method  is  shown  for  each  pairwise  comparison.  In  all 
cases,  the  T-test  showed  that  the  differences  in  the  mean 
error  were  significant  at  the  95%  confidence  level  or 
higher. 

Center  Errors 

Notice  in  Fig.  1  that  eddy  (c)  has  a  particularly  noncir¬ 
cular  shape.  For  all  methods  except  Janowitz,  this  caused 
a  spike  in  the  graph  even  if  it  was  not  the  largest  spike  for 
all  methods.  This  would  be  expected  for  detection  of  a 
noncircular  feature  using  a  circle-detection  method.  There 
seems  to  be  no  trend  for  better  detection  of  cold  over  warm 
or  vice  versa.  Of  the  method  tested,  Janowitz  is  the  most 
erratic.  The  existence  of  center  structure  and  large  missed 
edges  cause  problems.  Each  method  was  compared  to  the 
others  using  a  T-test  to  determine  statistical  significance 
of  the  result.  This  comparison  (Table  I)  shows 
Thomas_Chan  to  be  the  best  performing  of  the  methods, 
Albano  is  second,  and  the  simple  Circle  Fit  method  is 
third,  leaving  Hough,  Janowitz,  and  Albano  as  the  least 
effective. 

Radius  Errors 

Fig.  3  shows  the  results  of  each  method.  Albano  and 
Circle  Fit  show  trends  toward  defining  larger  than  analyst 
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Fig.  2.  The  six  methods  are  compared  for  center  position  with  the  analyst 
definition  of  the  18  eddies  (a)-(r). 


radius  parameters.  These  trends  do  not  appear  to  be  well 
defined  enough  to  apply  a  correction.  Hough  and  Janowitz 
seem  to  be  the  most  erratic.  Again,  the  methods  are  com¬ 
pared  using  a  T-test  (Table  II).  Again,  Thomas_Chan  per¬ 
formed  best,  with  Landau  second,  and  Albano  third. 
Hough,  Janowitz,  and  Circle_Fit  are  least  effective.  There 
would  be  no  problem  defining  center  position  and  radius 
by  two  different  methods,  but  the  results  show 
Thomas_Chan  to  perform  best  for  each.  The  eddies  used 
for  this  study  show  a  wide  range  of  characteristics— center 
structure,  missing  outer  structure,  noncircularity,  warm 
core,  and  cold  core.  The  data  set  is  small,  but  some  meth¬ 
ods  have  been  shown  to  fail  even  on  this  limited  data  set. 


Conclusions 

The  Thomas  Chan,  Landau,  and  Albano  methods, 
which  are  based  on  mathematical  minimization  of  some 
error  function,  tend  to  produce  the  best  circle  fits.  The 
nonmathematically  derived  methods  of  Hough,  Janowitz, 
and  Circle_Fit  tend  toward  poorer  results.  This  result  is 
not  necessarily  expected.  Certainly,  mathematically  op¬ 
timized  algorithms  will,  by  definition,  give  the  best  math¬ 
ematical  fii  to  the  detected  edge  fragments.  However,  it 


does  not  necessarily  follow  that  the  optimal  fit  to  the  frag¬ 
ments  will  in  fact  also  be  the  best  approximation  to  the 
judgment  of  a  human  analyst  in  marking  eddies  in  the  im¬ 
age.  Yet,  in  this  study,  the  mathematically  based  algo¬ 
rithms  also  most  closely  matched  the  subjective  human 
result. 

In  this  study,  the  Thomas_Chan  method  produced  the 
best  error  statistics  for  both  center  position  and  radius. 
However,  we  used  average  error  for  comparing  methods. 
Had  we  used  another  error  metric  (rms  error,  for  example) 
one  of  the  other  methods  might  have  been  the  best.  Also, 
circle-fit  results  will  depend  on  the  edge  detection  that 
precedes  the  circle- finding  algorithm.  Our  edge  detector 
is  typical,  but  if  another  edge  detector  had  been  used,  the 
ranking  of  the  methods  might  have  changed.  We  also  had 
only  18  eddies  in  our  data  set.  A  larger  or  different  data 
set  could  have  produced  slightly  different  results.  For 
these  reasons,  it  is  not  possible  to  definitively  select  a  best 
circle-fit  algorithm  based  on  the  results  of  this  study.  We 
can  conclude  that  the  mathematically  based  class  of  circle 
detectors  does  show  better  agreement  with  human  judg¬ 
ment,  even  when  the  machine  algorithm  is  presented  with 
only  fragmented  and  corrupted  pieces  of  the  circle,  than 
do  such  ad  hoc  algorithms  as  Hough,  Janowitz,  and 
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Fig.  3.  The  six  methods  are  compared  for  radius  with  the  analyst  definition 
for  the  18  eddies  (a)-(r). 


TABLE  I 

Best  Methods  Defined  by  T-Test  for  Center  Position 

Method  1  -  HOUGH  THOMASCHAN  LANDAU  ALBANO  JANOWITZ  CIRCLEFIT 

Method  2  i 

HOUGH  THOMASCHAN  LANDAU  ALBANO  HOUGH  CIRCLE_FIT 

THOMASCHAN  THOMASCHAN  THOMASCHAN  THOMASCHAN  THOMASCHAN  THOMASCHAN 

LANDAU  LANDAU  THOMASCHAN  ALBANO  LANDAU  CIRCLE_FIT 

ALBANO  ALBANO  THOMASCHAN  ALBANO  ALBANO  ALBANO 


JANOWITZ  HOUGH  THOMAS_CHAN  LANDAU  ALBANO  CIRCLE_F1T 

CIRCLE_FIT  CIRCLE_F1T  THOMASCHAN  CIRCLEFIT  ALBANO  CIRCLE_FIT 


TABLE  II 

Best  Methods  Defined  by  T-Test  for  Radius  Definition 

Method  1  -* 

Method  2  -l 

HOUGH 

THOMAS_CHAN  LANDAU  ALBANO 

JANOWITZ 

CIRCLE  FIT 

HOUGH  THOMAS  CHAN  LANDAU  ALBANO  JANOWITZ  CIRCLE  FIT 

THOMAS_CHAN  THOMAS  CHAN  THOMAS  CHAN  THOMAS  CHAN  THOMAS  CHAN  THOMAS  CHAN 

LANDAU  LANDAU  THOMAS_CHAN  LANDAU  LANDAU  LANDAU 

ALBANO  ALBANO  THOMASCHAN  LANDAU  ALBANO  ALBANO 

JANOWITZ  JANOWITZ  THOMAS  CHAN  LANDAU  ALBANO  CIRCLE  FIT 

CIRCLE_FIT  CIRCLE_FIT  THOMAS_CHAN  LANDAU  ALBANO  CIRCLE_F1T 
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Circ!e_Fit.  We  conclude  that  the  Hough  circle  detector 
that  v  incorporated  into  SAMAS  without  quantitative 
evaluation  was  not  the  best  choice. 

The  circle-fit  accuracy  in  this  study  was  approximately 
the  same  for  warm  and  cold  eddies.  This  result  was  a  sur¬ 
prise,  since  cold  eddies  offer  much  fainter  thermal  sig¬ 
natures  than  warm  eddies.  We  expected  that  the  weak  edge 
gradients  associated  with  cold  eddies  would  result  in  poor 
edge  detection,  which  would  in  turn  introduce  errors  into 
the  circle  fitting.  This  expectation  was  not  borne  out  in 
the  present  data. 

The  mean  position  error  is  1 1 .96  km,  and  the  mean  ra¬ 
dius  error  is  8.88  km  for  the  Thomas  Chan  method.  These 
error  values  provide  a  first  estimate  of  the  best  accuracy 
with  which  ocean  eddies  can  be  automatically  mapped  for 
IR  satellite  imagery. 
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